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EXPERIMENTAL DETERMINATION OF TIME CONSTANTS AND NUSSELT NUMBERS FOR 
BARE -WIRE THERMOCOUPLES IN HIGH-VELOCm AIR STREAMS AND ANALYTIC 
APPROXIMATION OF CONDUCTION AND RADIATION ERRORS 
By Marvin D. Scadron and Isidore Waxshawsky 

January 1952 


Page 20, equation (28): The numerical constant should be 0.0540 instead 

of 16.2. 

Page 20, equation (29): The numerical constant should be 0.00171 

instead of 0.246. 

Page 48, equation (C7a) : The exponent of T^ should be 3.815 instead of 

3.82. 

Page 48, equation (C7b) : The numerical constant should be 0.93x10"^® 

instead of 0.97X10-10 

Page 49, equation (C8b) : The numerical constant should be 4.05 instead 

of 4.20. 

Page 50, equation (Cllb) : The numerical constant should be 1.71X10“^ 

instead of 1.65X10-^ 

Page 69, figure 11(c): The ordinate values should be 10 times as large 

as the values indicated. 
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EXPERIMEHTAL I[EH?EBMINATION OF TIME CONSTANTS AND NUSSELT NUMBERS 
FOR BABE-WIRE THERMOCOUPLES IN HIGH-VELOCm AIR STREAMS 
AND ANALTTIC APPROXIMATION OF CONDUCTION 
AND RADIATION ERRORS 

By Marvin D. Scadron and Isidore Warshawslcy 


SUMMARY 


Experiments were conducted to determine time constants of several 
"bare-wlre thermocouples mounted In cross flow to, an air stream. Four 
pairs of common thermocouple materials as well as three different wire 
diameters were used. The prohes were tested at the exit of a room- 
temperature Jet exhausting Into a pressure controlled chamber. The Mach 
number range was 0.1 to 0.9; the range of Reynolds number based on wire 
diameter was 250 to 30,000. 

For the configurations tested, thermocouple material, wire diam- 
eter, gas pressure, and Mach number were shown to Influence the time 
constant. From these data, a correlation between Nusselt and Reynolds 
number was obtained for the range of Mach and Reynolds number stated. 

The correlation Is 




or 


Nu = (0.478±0.002)\/He* Pr°*^ 


with average deviations of a single observation of 6.9 and 7.4 percent, 
respectively. The value Re» Is a Reynolds number computed by evalu- 
ating gas viscosity and density at total temperature; Nu Is Nusselt 
number computed by evaluating gas thermal conductivity at total tempera- 
ture, and Pr Is the Prandtl number of the gas. 
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Graphs and nomographs are presented for the computation of approxi- 
mate radiation error, conduction error, and time constant for thermo- 
couples. 

Analytic determinations are also presented of the effects of dis- 
similarities in vire material and wire diameter on the steady- state tem- 
perature distribution in the exposed portion of the wire, and of the 
effect on time constant of conduction along the thermocouple wire to the 
supports . 


HJTEOnJCTION 

Perfomance evaluation and control of jet engines, as well as the 
fundamental study of the related combustion phenomena, depend consider- 
ably on a knowledge of the steady and variable temperatures of the gas 
within the engine. At present, these temperatures are most commonly 
measured with thermocouples. The measurement accuracy is generally 
dependent on the temperature level, the gas velocity, the tenperature of 
the surrounding walls, and the thennocouple construction. The design of 
a thermocouple for any particular application represents a compromise 
among the contradictory factors of accuracy, life, ruggedness, and 
rapidity of response to changes as they are influenced by conduction and 
radiation losses, partial adiabatic recovery, erosion, size, and con- 
vective heat -transfer rate with the moving gas. 

Whether a thermocouple is designed to respond rapidly to tempera- 
ture changes or to measiire accurately the steady- state temperature of 
the gas, the controlling factors are the same, although the orders of 
importance may be different. Although this report is concerned prin- 
cipally with speed of, response, consideration of the related steady- 
state acctaracy is included. 

Limited data have been available for estimating time constants and 
these data are for certain thermocouple designs and for a smiall range of 
operating conditions. Data presented in reference 1 are time constants 
of becre-wire and radiation- shielded thermocouples at Eeynolds numbers 
(based on wire diameter) from about 15 to 900 and Mach numbers from 0.05 
to 0.14. The work was performed in exhaust gases at 1000° F, at gas 
velocities up to 250 feet per second, and with a step c h a n ge of approxi- 
mately 700° F. For a given thermocouple, the time constant was found to 
vary with the 0.5 power of the Beynolds number. Additional response- 
rate data on very fine wires for a Eeynolds number range of 3 to 340 and 
a Mach number range of 0.02 to 0.1 are presented in reference 2. As 
shown in reference 2 and in other references, the study of response 
rates can bo tantamount to the measurement of Nusselt number. The rela- 
tion between Eeynolds number and Husselt numOber, and therefore response 
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rate, for cylinders in cross flow is represented in a compilation of 
data in reference 3 for a Eeynolds number range of 0.1 to 250,000 and 
Mach numbers up to 0.1. As gas velocity approaches the sonic value, it 
is conceivable that compressibility effects, as represented by the Mach 
numiber, may influence the Nusselt number. This influence is shown by 
the data presented in reference 4. These data were obtained in the 
Eeynolds number range of 80 to 550 and Ifech number range of 0 . 4 to 2.2. 

The present work was undertaken to evaluate time constants of some 
representative types of thermocouple probe that are applicable for use 
in Jet-engine operation. Experimental determinations were made of the 
time constants of bare-wire thermocouples in cross flow to an air stream 
for Eeynolds numbers between 250 and 30,000 and for Mach numbers between 
0.1 and 0.9. Nusselt numbers were computed from the experimental meas- 
urements. The ranges of Eeynolds and Mach number reported by refer- 
ences 1 to 4, inclusive, as well as those presented herein, are shown in 
figure 1. 

The associated theoretical analysis serves to evaluate the possible 
effects of conduction on the time constant and to provide an estimate of 
the conduction and radiation errors that can be expected in steady- state 
temperature indication. For the range of conditions covered, the data 
aid in establishing the relation among Nusselt, Reynolds, and Mach 
numbers to an accuracy adequate for estimating conduction and radiation 
errors and thermocouple time constants. 

The work reported herein is part of a program of research in hlgh- 
temperatoire measurements being conducted at the NACA Lewis laboratory. 


THEORT 

Method of theoretical treatment . - The theoretical treatment of 
this section will first define the concept of "time constant" as it 
appears in the most elementary case of convective heat transfer. Then, 
in the treatment of the more general case of combined convective, con- 
ductive, and radiative heat transfer from a single wire, the manner in 
which the time constant is affected by the presence of conduction and 
radiation will be shown. This treatment will lead to theoretical form- 
ulas for the approximate steady-state conduction and radiation errors. 
The interrelation between the thermocouple response rate and the factors 
affecting the steady-state errors will be demonstrated in the section 
Basic example. Next, there will be considered several more complex 
situations of practical interest: the single wire replaced by a thermo- 

couple, a thermocouple with intermediate supporting wires, a single wire 
subjected to a step change from initial conditions different from those 
considered under Basic example, and a single wire heated electrically. 
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Finally, It will te shown how the introduction of en^jirically determined 
constants into the theoretical formulas may be used for the numerical 
computation of the approximate values of radiation error, time constant, 
and conduction error. The symbols used throughout this report and the 
detailed mathematical procedures that lead to the results presented in 
this section are given in appendixes A, and B and C, respectively. 

The heat- transfer relation for a bare-wire thermocouple used to 
measure the temperature of an air stream will be developed under the 
following assumptions : 

(a) The difference between thermocouple temperature and air- stream 
temperatiare is small compared with the absolute temperature level. 

(b) The thermocouple dimensions are small compared with the dimen- 
sions of the duct or other solid enclosure. 

(c) Eadlal temperature gradients in the wire may be neglected, so 
that the temperature at any cross section is constant throughout that 
cross section. 

(d) The thermocouple is constructed in one of the forms shown in 
figure 2, with the thermocouple junction midway between the supports, 
and with the support temperature remaining constant or else changing 
very slowly as compared with the thermocouple- Junction temperature. 

Elementary case, excluding conduction and radiation . - In the 
absence of conductive and radiative heat transfer, the balance between 
rate of accumulation of heat in an element of len^h of wire and rate of 
convective heat transfer to the same element of length would be stressed 
by a simple total differential equation 


V. ■^-dJ - “<’^8 - 

where 

A surface area 
c^ specific heat of wire 
h convective heat- transfer coefficient 
T effective gas temperature 
T^ wire temperature 
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t time 


V volxime 

wire density 

Equation (la) may be written 


T 



+ Tv 


‘•g 


(lb) 


where t is the time constant of the wire. Characteristic of equa- 
tion (lb) are the facts that 

(a) In response to a sudden "step change" in gas temperature from 

an initial value (at which time will be taken equal to T^^^) 

to a new constant value Tg "tli® wire response will be given by 

■'w ' ■'s.2 * 

(b) In response to a sinusoidal fluctuation in gas temperatxire, of 
angular frequency ui, represented by 


T = T +T sin ((ot) 
g g,av ^ ^ ' * 


and initiated at time t = 0 (at which time will be taken as equal 

to Tg g^^) , the wire temperature will be given by 


T 


w 


g,av 


+ sin 




T QQT 

1-KU^ 


,-t/T 


(3) 


General case . - In the presence of conductive and radiative heat 
transfer, equation (l) must be replaced by a partial differential equa- 
tion involving z, distance along the wire, as a second Independent 
variable. The solutions for response to a step change or to a sinus- 
oidal fluctuation are correspondingly more complex. Examination of 
these equations and comparison with equations (2) and (3) indicate, 
however, that, under the assumptions initially stated, it is still poss 
ible to have a characteristic time constant t for the thermocouple 
(although it may no longer be numerically equal to that appearing in 
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equation (ib)) and the behavior of the thermocouple junction can still 
be described by equations of the form of equations (2) and (3). 

In order to indicate the form of the time constant in the more 
general case, as well as to evaluate the errors in steady-state ten 5 )era- 
ture indication that are caused by conduction and by radiation, the gen- 
eral heat-transfer equation will be derived. 

For an element of length Ax of one of the wires, the rate of 
storage of heat q_ is equal to the sum of the rates of heat transfer 
into the wire element by the processes of conduction, convection, and 
radiation. These rates of heat transfer will be denoted by q}^, q,,, 

and qj., respectively, so that 


qp Ax = (qj. + + qj.) Ax 


U) 


The following additional symbols will be used: 

I 

D wire diameter 


g 


thermal conductivity of gas 


w 


thermal conductivity of wire 


L length of wire between supports 


Uu Hus Belt number 

Tj equivalent duct temperature, defined as absolute teii 5 )erature of a 
black-body enclosure (of included solid angle 4n) whose radia- 
tion to thermocouple, through gas, will be same as that of 
actual duct 


«g 

€w 

0 


absolute wire temperature at thermocouple Junction 

effective absorptivity of gas for black-body radiation at 
temperature 

effective emisslvity of gas 
emlssivity of wire 
Stefan-Boltzmann constant 
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The rate of storage of heat in unit length of the wire is 


% 


= Pw 


c 

^ ^ 4 


( 5 ) 


The net rate 
unit length of wire 


of heat transfer by conduction through the ends of 
is 




S^T^ kT? 


(6) 


The rate of heat transfer by convection into the surface of the 
wire per unit length of wire is 

Nu kff , , , . 

■’c - “7^ fV'-’ 

From the surface of unit length of wire, the radiant power emitted 

^Cglg^xD, and the 

power received from the duct by transmission throu^ the gas is 
cre■^,(l-ag ^)T< 3 i^jcD, Hence, the net rate of heat transfer by radiation to 
the surface of unit length of wire is 

9r - os, - ’’w^] 

However, if approximation Tg^-T^^ ® 4Tg^(Tg-T^) 

may be introduced, permitting equation (8) to be written 

An analysis of the shape factors that enter into the determination of 
T^ is presented in reference 5. 


is the power received from the gas is o€ 
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Substitution of equations (5), (6), (7), and (9) into equation (4) 

yields 


T 



1 


— + If - T. 


-w 


( 10 ) 


where t, sre defined by 

4 Nu kg 


(11a) 


h 



Hu kg 

4 Nu k„ 
d 2 k^ 


T 


XiSepTij 


(lib) 


(lie) 


(lid) 


3 = 


- PlS. 

(l+43^€^/Tg) 


(lie) 


^ (l+43i€^/Tg) 


(Ilf) 


Tf - Tg + 3 




n-..g 


(llg) 


For purposes of future computation, it is also convenient to 
introduce the thermal diffusivity of the wire 


K = 


^w 


P c 

w 


(llh) 
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so tliat equation (llf) may also "be written 


* i/k t 


(Hi) 


The physical significance of the quantities appearing in equa- 
tions (11) hecomes apparent upon examination of the solutions of equa- 
tion (10) for some simple boundary and. initial conditions, as treated in 
appendix B. 

Under conditions of a steady gas temperature, the corresponding 
steady state of the wire =0) is one in which the balance of 

heat-transfer rates by conduction, convection, and radiation produces a 
wire temperature that is different from the gas temperature Tg 

and is given by 

= Tg + p|jl-ag^4)^^^ - (l-€g)^ (Tb-Tf)iJ/ (12a) 

where * >^(ti,x), and T^ is the temperature at the supports, x * 0 
and X a: L. 

The quantity T includes the effect of radiant heat transfer and 
represents the temperature the wire would asstnne if there were no con- 
ductive heat transfer. The "radiation error" is therefore given by 

An alternative form for the radiation error is 

where 


(l+4Pi€^€g/T;,) 


(13a) 


(13b) 


(13c) 
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and 


Pi 


gP 

m kg 


(13d) 


This alternative form is useful vhen the indicated wire tenqperature 
is known, and it is desired to determine the radiation correction to he 
applied to the reading in order to obtain Tg. The use of the alterna- 
tive quantity "p permits expression of equation (iHa) in the form 


T 


g 






(1-Cg) 


(T^-T^) -if 

1 -, 1 , 


(I2h) 


The intermediate quantities Pi and are introduced merely for con- 

venience in writing the equations. The effect of conductive heat trans- 
fer is to make the steady- state wire temperature different from ^T^. 

The magnitude of this difference and hence the "conductive error" is a 
function of q and is directly proportional to (Tj.-T-jj). An additional 

effect of conductive heat transfer is to produce a nonuniform tempera- 
ture distribution along the length of the wire, ranging from a value T^ 
at the supports to a value T^ ^ at the midpoint of the wire (thermo- 
couple junction) . ’ 

If both thermocouple wires have substantially the same diameter and 
thermal conductivity, they also have substantially the same q. The 
steady-state conduction error at the midpoint, where x » L/2, is then 


(V^f ) 'I'm 


(14) 


where, as shown in appendix B, section II, 


= sech (qL/2) 


(15) 


If the gas temperature is changed to a new value, the wire 
approaches a new steady- state value, involving a new value of T^, The 

approach to the new value is at a finite rate, determined by the quan- 
tities T and q. In the absence of appreciable conductive and radia- 
tive heat transfer, the wire has a time constant in the absence 

of conduction only, the time constant is t; in the presence of a 
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moderate amount of conduction, the thermocouple Junction acts approxi- 
mately as though it had a time constant 


T 


(16) 


Basic example , - An illustrative example is the simple case vhefe . 
both thermocouple wires have substantially the same diameter and thermal 
conductivity. If (a) the gas temperature is suddenly changed from a 
value T T to a new steady-state value T (step change), (b) the 

thermocouple was initially in equillhrium with the gas, and (c) T , > 

S/ i 

Tg^2 "^dj response of the thermocouple junction is shown in 

figure 3(a), where the yariatlon of with time is Indicated. The 

initial and final distributions of temperature along the wire are shown 
in figure 3(b), where the variation of with is indicated at 

times t = 0 and (t/r) *. The derivation of the equation for steady- 

state temperature distribution along the wire is presented in section II 
of appendix B. The effect of parameter qL upon the distribution is 
shown in figure 4. The equation for time response of the thermocouple 
Junction is derived in section VT of appendix B, resulting in equations 
(B39) and (B40). If the wire is subjected to a sinusoidal variation of 
temperature rather than a step change, the time response of the thermo- 
couple Junction, derived in section VII of appendix B, is given approxi- 
mately by equation (B53) . 

Comparison of equations (B40) and^BSS) with equations (2) and (3), 
respectively, leads to the conclusion that the effect of conduction is 
to attenuate the temperature change by a factor of (l-\|(m), to multiply 
the time constant t by the same factor, and to add a steady conduction 
error equal to times the difference between the support temperature 

and the average impressed temperatiure 

Other examples . - In' addition to the simple example Just discussed, 
four basic alternative situations are of interest. The first of these 
is the case when the two thermocouple wires have radically different 
values of the quantity (k^riD^). The expressions for conduction error 
and effective time constant become more complex, but an equivalent 
quantity (n'L) can be defined in terms of the quantities (tj^L + Igb), 

(Ha^ ~ tIbL), and (k^TiD2)^/(]c^TiD^)g, where the subscripts A and B 

pertain to the two thermocouple wire materials. An equivalent quantity 
'If'm can be used in the expression for steady-state and transient res- 
ponse discussed previously. The steady-state conduction error for this 
situation is discussed in greater detail in section III of appendix B. 
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Tli3 sscond fil'tomfl'tiv© si’tus'tioiL of in'torss't Is oo© in wiilcli "tii© 
thermocouple Junction Is made between two wires which are attached to 
the supports throu^ intermediate wires of diameters different from 
those forming the Junction (fig. 2(c)). Here again equivalent quanti- 
ties (q"L) and can he defined in terms of the tjL values and 

the ky ratios. The steady-state conduction error for the case 

where both thermocouple materials have substantially the same thermal 
conductivity is treated in detail in section 17 of appendix B. Combina- 
tion with the case discussed in the preceding paragraph Is possible in 
order to include the effect of an appreciable difference between the 
thermal conductivities of the two wires. 

The third alternative situation is one in which a step change to a 
constant value of gas temperature Tg^g applied to a thermocouple 
that is initially at support temperatiire T, rather than being initially 
at equilibrium with the gas stream. This pfocedure is descriptive of the 
lag- testing technique described in reference 1. The thermocouple 
is initially shielded from the hot gases by a tube through which cooler 
air is blown. The tube is then suddenly removed, exposing the unit to 
the hot gases. The response is shown graphically in figure 5. Mathe- 
matically, it involves replacement of the multiplier (Tf^q-Tf^s) in 
equation (B40) by the multiplier (T^j-Tf^g)* 

The fourth alternative situation is represented by the esperiments 
described in this report in which the heating effect of an electric 
current was used to raise the temperature of the wire artificially and 
the step change was produced when the current was broken. This situa- 
tion involves the addition of a term qeAx = W Ax to the right side of 
equation (4), where W is the power dissipated per unit length of wire. 
The effect of this term is merely to add the quantity W/rt Hu kg to 

the right side of equation (llg) for Tj. The effective time constant 
is unaffected. The response is shown graphically in figure 6. 

Computation of time constant, conduction error, and radiation 
©j^orT —The Nusselt number Nu enters into formulas (11a) through (llf) 
for the constants t. , 3q, and that are required for the com- 

putation of radiation error, conduction error, and effective time con- 
stant. It is conventional in heat-transfer work to express this number 
in terms of the fluid properties and dynamic conditions by means of the 
relation 


Nu = F(Be, Pr) 


(17a) 
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where Ee and Pr are the Eeynolds number and Prandtl number, respect- 
ively. Because of the large range of aerodynamic conditions covered in 
these tests, the conditions being representative of the possible uses of 
the thermocouples, it was found appropriate to modify equation (17a) to 
include the Mach number M as a separate independent variable and to 
use a simple product of powers of the Eeynolds and Prandtl numbers as 
part of the function F, thus yielding , ' 


Nu = constantxEe® Pr^ i(M) 


(17b) 


For the purpose of actual computation, it is desirable, so far as 
is convenient, to express the values of 3^, and in terms of 

quantities that are actually measured. Such quantities are the total 
temperature and the static and total pressures. For greater conven- 
ience, the Mach number M is also Introduced as an Intermediate vari- 
able, its value in terms of static- to total-pressure ratio being read- 
ily available in tables. Then, utilization of the additional facts that 

(a) both the viscosity and the thermal conductivity of air can be ex- 
pressed as proportional to a power of the temperature (appendix C) and 

(b) the Prandtl number for air probably enters as the 0.3 power and the 
value of Pr®*3 may be treated as remaining constant (appendix C) shows 
that the Eeynolds number and, consequently, the Nusselt number can be 
expressed as products of powers of the Miach number, static pressure, and 
static temperature. In appendix C, the eapressions for Eeynolds number 
and for thermal conductivity of air are shown to be 

-1.19 I— 

Ee = const ant XDMpTg ■ VT (18) 

kg * constantxT^*^® (19) 


After these expressions are inserted into equation (l7b) and the static 
temperature is used in the equation for thermal conductirity^ equations 
(11a), (lib), and (lie) may be "written 


T"! * 


2-a 

p c D 

V 

paj^ar O- 78-1.193 ^a/z 
s 


constantx- 


(20a) 
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p®M®T f(M) 

= constantx ^ — 

ky-= 


p, = constantx 

p^M® f(M) 


(21a) 


(22a) 


where the Frandtl number has been included as a part of the constants. 

If gas density, viscosity, and thermal conductivity are evaluated 
at total rather than at static taaperature , an alternative Reynolds 
number Re* is defined 


Re* « constantXDMpT^ ^ xj 

(appendix C) and alternative forms of equations (I7b), (20a), (21a), and 
(22a) are 


Nu » constantx(Re*)® Pr^ f*(M) 


l/2 


■ constantx- 


pVr 7®/^ f*(M) 

X 


(17c) 

(20b) 


constantx- 


j »^^ 0 . 78 - 1 . 19a ^./2 




i/2 


3. 22+1. 19a, 


aD T. 


P » constantx- 




i/2 


p®M® 7^^^ f*(M) 


(21b) 


(22b) 
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APPARATUS AUU PROCEDURE 

•PT-l ph-veloclty air source . - The tests were performed at the exit of 
a 2^- inch-diameter jet that exhausted into a pressure controlled re- 
ceiver (fig.’ 7). Performance characteristics of the Jet are described 
in reference 6. At the test section, a Mach number of 0.1 to 0.9 and a 
free-stream density of 0.04 to 0.10 pound mass per cubic foot were avail- 
able. By individual control of the Jet inlet and receiver pressures, 
any combination of Mach number and density for the ranges listed was 
obtainable. Each run was made at constant test-section density and 
varying Mach number. The total temperature of the air was between 540° 
and 545° R and was constant to within 1° during any one reading of a 
run. The total pressure and total temperatxire in the plenum chamber, 
and the static pressure at the test section, with the probe inserted 
into the air stream, were used to calculate Mach number and free-stream 
density. 

Test probes . - Six thermocouple probes were tested; their dimen- 
sions and configurations are shown in figure 8. The foTir pairs of wire 
materials used were_ chromel-constantan, chromel-alumel, Iron-constantan, 
and platinm plus 13-percent rhodium - platinum. Probes A, B, and C 
were of the same wire material but of increasing wire diameter; probes 
B, D, E, and P were of approximately the same size and configuration but 
of different wire materials. A seventh thermocouple was also tested for 
purposes of comparison with probe B. This thermocouple. Identified by 
the letter G, was a straight chromel-constantan bare wire stretched 
across the exhaust of the Jet. 

Time- constant determinations . - The time-constant determinations 
were made without moving the probe or shielding it from the air stream. 

In order to simulate a higher initial gas ten^jerature, an electric 
current was passed through the thermocouple wire to heat the Junction to 
about 100° R above room temperatiore. A double-pole, double- throw snap 
action toggle switch was then used to disconnect the thermocouple from 
the source of heating current and to connect it to a recording oscillo- 
graph. As Indicated under TBEQRY, this technique served to simulate a 
step change from one gas temperature to another. It also served to 
minimize any change in the nature of the ten^ierature distribution along 
the wire. This distribution itself was, for all wires tested, such that 
the contribution of the term 'Hi/ to the effective time constant was 
negligible. For this reason, ai^d because of the technique used for 
measuring the time constant as described later, it was also possible to 
neglect without serious error two deviations from the assumptions of the 
simple theory. These deviations were: first, the base of the thermo- 

couple loop, constituting the "supports," actually did change tempera- 
ture moderately; and second, for the iron-constantan wires an appreci- 
able difference exists between the electric power dissipation in the two 
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halves of the thermocouple loop and, consequently, an asTnimetry In tem- 
peratiire distribution and a difference in the values of occtot. For 

exaizq)le, in the most serious case, that of probe D at M = 0.1 and 
p » 0.04 pound mass per cubic foot, the value of qL is 6.6 for the 
iron vlre and 11.0 for the constantan wire, corresponding to a correc- 
tion of 3 percent if the support temperature is assumed equal to the 
initial wire temperature. The ratio of tjL values for the two halves 
of the wire was 1.7:1 and the ratio of increase in values due to 

the heating current was 5:1. 

As the thermocouple was cooled, ita output was recorded by a criti- 
cally damped moving-coil oscillo^aph element hav ing a natural period of 
0.006 second and hence an effective time constant of 0.002 second. The 
record was obtained upon photographic paper that was marked off in 0.01- 
second intervals. The paper speed and oscillograph-element sensitivity 
were constant for all practical purposes. A typical record is shown in 
figure 9(a). Figure 9{b), a graph of the same data on semilogarithmlc 
paper, shows that in the region following the "Start of record" the 
relation is sufficiently linear that a single measurement, directly on 
the oscillograph record, of the time required for the recorded tempera- 
ture to have covered 63.2 percent of the tenq>erature interval from start 
to asynqptote is sufficient to establish the time constant for each of 
the runs. The time constants were read to a probable error of 0.005 
second. For the temperature range encountered in these experiments, it 
was reasonable to assume that the wire calibration (emf/deg temperature 
change) remained constant. Each of the measured values of time con- 
stant was used as obtained from the records since the maximum conduction 
error by equation (16) was never more than 3 percent and usually about 
1/2 to 1 percent. 


EESULTS AND DISCUSSION 

T^we constants and Nusselt numbers . - The data obtained in the 
experiments are shown in figure 10 as a plot of the experimentally meas- 
urecTtime^ constants against the Eeynolds number Be based on static 
pressure and temperature. The conduction correction factor is 

negligible for almost all the data, being no more than 3 percent in the 
extreme case. The radiation correction is also negligible. Conse- 
quently, the experimentally measured time constants are assumed to be 
the same as the quantity of equation (lla). 

It will be noted that points in the lower Mach number region fall 
on 3-ines having a constant slope of approximately 0.5 and that points in 
the higher Mach number region show s^tematic deviations from these 
lines of constant slope. This effect is more strongly Indicated by 
plots of time constant as a function of Mach number, as shown in 
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figure 11. On the average, these figures indicate that the deviations 
begin at approximately M = 0,5. 

An alternative presentation of the data is a plot of the time con- 
stant against the Eeynolds number He* based on an, evaluation of gas 
density and viscosity at total temperattire rather than at static temper- 
ature. Such plots (fig. 12) show appreciably smaller deviations at 
higher Mach numbers than were in evidence in figure 10 and ars equivalent 
to assuming for the functions f(M) and f*(M) in equations (17b) and 
(17c) the forms 

•v-1 o “I'SOa 

f(M) = (1 + 21— m 2) (23a) 

f*(M) = 1 (23b) 


These expressions follow from the fact that 

-Re* . Ee\(l + Zll M^) (24) 

2 


(See appendix C.) If Ee* replaces Re, it is reasonable to assume 
that the total temperature, rather than the static temperature, should 
be used in evaluating the thermal conductivity. Making this assump- 
tion and utilizing the values of mean Pv^w tables I or II and 

of thermal conductivity of air at the operating total temperature per- 
mit computation of the Nusselt number by equation (lla) . A plot of 
this computed Nusselt number against the Reynolds number Re* is 
shown in figure 13 for all the 176 data points taken in these tests. 
The least-squares solution for the best straight line throiigh these 
points, assuming that all Nusselt number measurements have equal 
probable percentage error, is 


Nu = (0.385*3.016) (Rer») 


(0.515±0.005) 


(25a) 


in the Eeynolds number range 250 < Re* < 30,000. The average deviation 
of any point from the line is 6.9 percent. If the exponent of Ee*- is 

fixed at the value i, the least-squares solution for the multiplying 

2 

constant leads to 


Nu = (0.431*0.002) ^Ee» 


(26a) 
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vitli an average deviation of any point of 7.4 percent. Fommla (26a) 
will be used in future discussion because of its greater convenience and 
negligibly greater average deviation. 

If tbe Husselt number - Beynolds number relation is to be applied 
to gases other than air,, the Prahdtl number Pr may be inserted explic- 
itly in equations (25a) and (26a) by writing them 

Nu = <0.427±0.018)(Ee*)^°‘^-^^*°°^^ Pr°*^ (25b) 

Nu = (0.478db0.002) Pr°*^ (26b) 


where the exponent b in equation (l7b) has been set equal to 0.3 in 
accordance with the conclusions reported in reference 3, chapter VIII. 
Formulas (25a) and (26a) are Independent of any assunqitions regarding 
the value of Prandtl number; formulas (25b) and (26b) are based on a 
value of Pr » 0.71 for air at 540° E. 

/ 

Some of the sources of error, in estimated order of Importance, are 
lack of roundnesB and smoothness of the wire. and Junction, error in 
effective wire-diameter measurement, deviation from simple theory be- 
cause of changing support temperature. Inaccuracy in measurement of time 
constant, uncertainty in knowledge of p^c.^, use of the arithmetically 

averaged PyC.^ for the two thermocouple materials, and omission of the 
conduction correction. The accuracy of the correlation represented by 
equations (25) and (26) is dependent on the accuracy of knowledge of gas 
viscosity and thermal conductivity at the operating total temperatxire 
(540° B); because all tests were made at constant total temperature, the 
accuracy of the correlation and therefore formulas (25) and (26), are 
independent of assumptions made concerning the variation of Prandtl 
number, gas viscosity, and thermal conductivity with temperatiire. 

A few experiments performed with probe C at a gas temperature of 
1500° E showed that the formulas presented remain applicable. The data 
are not presented herein. 

The compilation of data presented in chapter VIII of reference 3 
gives, for the range M < 0.1, 1000 < Be < 50,000, 

0.60 


Nu - 0.26 (Be) 
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Data (at hi^er gas temperatures) from reference 1 give the value of 0.5 
for the exponent of fieynolds number for the range M < 0.14, 20 < Ee < 
1000. Data from reference 4, at M = 0.38, give 

Nu = 0.47 ^|ie 


The data in references 3 and 4 were based on an evaluation of gas prop- 
erties at a temperature substantially equal to the arithmetic average 
of wire temperature and total gas temperature. For purposes of compari- 
son, the lines represented by formulas (25a) and (26a) are shown in 
figure 13 together with the curve given in reference 3. 

Effective gas temperature T„. - Most satisfactory correlation 

— -- - - O 

between time constant and Reynolds number is achieved by evaluating gas 
density and viscosity at total temperature. This technique can be ratlon- 
ed-ized by considering that the wire temperature Ty substantially repre- 
sents the average ten^ierature of the gas immediately surrounding the wire 
and that this layer of gas offers the principal resistance to the flow of 
heat. A graph of recovery temperature ratio T^/t^ of a cylinder in 
cross flow as a function of Mach number is shown in figure 14. The data 
are the composite results of a large number of experiments performed at 
this laboratory during the past 5 years in connection with other research 
projects. These measurements were made so that conduction and radiation 
errors were minimized and therefore Tg would be approximately equal to 
Tw for these data. These data agree quite well with values reported 
from time to time in the scientific literature. The cross-hatched 
section on the graph indicates the average spread of all the available 
data, which include the effect of different free-stream air densities. 

The ratio of T^/T^ is also shown. It is evident that the recovery 

temperature T^ is much closer to the total temperature Tt than to 
the static temperature Tg. Reynolds numbers computed at M = 0.9 based 
on gas properties evaluated at T-t, Ty, and Tg woxild be in the ratio 
of 1.00:1.00:1.29, respectively. 

The finding that a batter correlation of heat-transfer data is 
often achieved when gas properties are evaluated at surface temperature 
was first reported in reference 7 for low-velocity flow and has recently 
been verified at higher velocities in reference 8. The experimenters in 
references 7 and 0 were concerned with flow in pipes, where viscosity is 
the predominant factor, rather than with flow around cylinders, where 
potential flow effects are as important as viscous effects. Conse- 
quently, the physical explanation of the finding that the use of total 
temperature provided a better correlation than the use of static temper- 
ature in the case of flow around cylinders is not necessarily of the 
same nature as the explanation of the heat- transfer phenomena in pipes. 
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HDMERICAL CCMPUTATION OF TIME CCWSTAHT, 

RADIATION ERROR, AND CONDUCTION ERROR 

Tlie Nusselt number, once having been dotexTiilned from experimental 
measurements, may then be inserted into the basic theoretical formulas 
for time constant, conduction error, and radiation error. Since, as 
shovn in appendix C, air viscosity and thermal conductivity may be ex- 
pressed as functions of absolute teniperature, it follows that Reynolds 
number, Nusselt number, radiation error, time constant, and conduction 
error, may also be expressed in terms of static pressxire, total pres- 
sure, and total temperature, all of which are experimentally measTirable 
quantities, Tbe computation procedure is facilitated by introducing 
Ifach number, the value of which is readily obtained from tables if the 
static and total pressures are known. 

The insertion of equation (26a) into equations (20b), (21b), and 
(22b) leads to the equations 


4.05 (x X 

j - 

pO.50 yO.SO T^O.ia 


seconds 


(27) 


^1 = 


16.2 a D 


0.50 


25 


0.50 j^.SO 


°R 


^1 


0.246 p°-®° 


\0.25 


where the quantities are in the following units: 


Inches 


-2 


(28) 


(29) 


Btu/(ft2 Or) 

D (in.) 

p (atm) 
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Btu/(ft °R sec) 
a 0.173X10"® Btu/(hr ft^ 

These fontrulas show that time constant and conduction error are only 
slightly affected by the magnitude of the gas temperature; hence the 
validity of these formulas is not strongly dependent on any assun 5 )tions 
made concerning dependence upon ten^ierature of the gas viscosity and 
thermal conductivity. The time constant con 5 >uted from equation (27) has 
been plotted in figure 15 against the experimentally measured time 
constants. As 1s to be expected from the fact that figure 15 is sub- 
stantially an alternative presentation of the data shown in figure 13, 
the average deviation of the points shown in figure 15 la also 6.9 
percent. The detailed derivations of numerical formulas such as equa- 
tions (27), (28), and (S9) are presented in appendix C and lead to the 
construction of graphs and nomograms (figs. 16 to 19) that materially 
facilitate the computations of time constant, radiation error, and con- 
duction error for wire diameters between 0.001 and 0.1 inch. The exact 
computation procedure follows: 

Given: 

(1) Wire-material constants: €^, p^, c^, or K 

(2) Probe constants: D, L, L*, Tjj, T^ 

(3) Constant of surroundings: T^ 

(4) Gas and aerodynamic constants; M, p, a.^ 

Procedure ; 

(1) Find from figure 16(a). 

(2) Compute p = 3^ €^/(l+4B3^€^£g/T„) 

(3) Conpute radiation error => "^(l-a^ ^g^] 

the nomograph of figure 16(b) . 

(4) Find Tp^ from figure 17. 

(5) Compute = Tp^ (p^^^ c^^^ + p^^g c^,B^/®Ppt°Pt 

using data in figure 17. 

(6) Compute r =• T-j_/( 1+401 e^^g/Tv) 
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This value is the time constant in the absence of conduc- 
tion effects. 

(7a) Compute the effective (ilL) =» L / t K using the harmonically 
averaged value of k from table II. 

(7b) Alternatively, find (t]^L)j>^^ from figure 10 and compute, 
for each thermocouple wire, 

and 

(iL)^2 ^ (ti^L)„^ (l+4PiC^€g/T^) 

where (kp^/k^) is listed in table I, Then find the 
effective (ttL) for the thermocouple wire pair by use of 
the equations 

V”b ■ 

"a,b = [(V’‘b>+^] 

and of tables II and III. 

(8) Find from figure 19. 

(9) Compute steady-state conduction error = (T^ - T^) 

(10) Canpute corrected time constant * t (l - \j/ ) 

in 

(11) Compute corrected gas temperature Tg = - (steady-state 

radiation error) - (steady-state conduction error). 

SUMMAHT OF RESULTS 

The time constants of hare-wire thermocouples mounted in cross flow 
to an air stream were shown to depend principally on thermocouple 'wire 
material, wire diameter, gas pressure, and Mach number, and only 
subtly on gas temperature. These measurements also served to provide 
a correlation among Mach number M, Nusselt number Nu, and Reynolds 
number Be* for the ranges 0.1 < M < 0.9 and 250 < Be* < 30,000, 
yielding 



HACA TN 2599 


I 


23 


^ ,(0.515^0.005) 

Nu = (0.385i0.016)(Ee»)' 

with a 6,9 percent average deviation of a single observation, or 

Im =■ (0.433±0.002) 

with a 7.4 percent average deviation of a single observation. The value 
Eo* is a Beynolds number computed by evaluating gas viscos.ity and dens- 
ity at total temperature. The value of Eu was computed from the 
experimental measurement of time constant and the evaluation of gas 
thermal conductivity at total temperature. Prandtl number was consid- 
ered constant for these tests and is therefore Incorporated in the con- 
stants. The Prandtl number Pr may be Inserted explicitly in the pre- 
ceding equations by writing them 


/ ^ W ,(0.515±0.005) „ 0.3 

Eu » (0.427±0.018)(Ee»)' Pr 


Eu - (0.478±0.002) ^e* Pr 


0.3 


It is shown that the evaluation of Eusselt number provides a means 
for determining approximate steady-state radiation and conduction errors 
of a bare-wlre thermocouple in a hl^-temperatxire gas stream. Graphs 
and nomograms are presented for the computation of time constant^ 
approximate radiation error^ and approximate conduction error for wire 
diameters between 0.001 and 0.1 inch for commonly used pairs of thermo- 
couple materials. 

Analytic determinations are presented of the effects of dissimilar- 
ities in wire material or wire diameter on steady-state temperature dis- 
tribution in a thermocouple^ and the effect on the time constant of con- 
duction along the thermocouple wire to the support. 


Lewis Plight Propulsion Laboratory 

Eatlonal Advisory Committee for Aeronautics 
Cleveland, Ohio, August 3, 1951 
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APPENDIX A 
SIMBOLS 

The following symbols are used In this report; 
surface area 
specific heat 
wire diameter 
heat-transfer coefficient 
thermal conductivity- 
length of wire between supports 

length of central wire when intermediate supporting wires are 
used 

free- stream Mach number 
k D^ Tj 

Nusselt number 
Prandtl niimber 
static pressure 

rate of heat transfer by convection, per unit length 

rate of heat transfer by electric power dissipation, per unit 
length 

rate of heat transfer by conduction, per unit length 

rate of heat storage, per unit length 

rate of heat transfer by radiation, per unit length 

Beynolds number with gas properties computed at static tempera 
ture 
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Be*" Eeynolds number witb. gas properties con^juted at total temperature 

r specific gas constant for air, 1716 f b^/(sec^) (°R) 

T temperature 

T amplitude of sinusoidal gas -temperature fluctuation 

equivalent duct temperature 

T final wire temperature in absence of conduction 

f 

t time 

U linear free- stream velocity 

u intermediate temperature variable 

V volume 

V Intermediate temperature variable 

W electrical power per unit length 

V Intermediate temperature variable 

X distance along wire 

T,y,z intermediate temperatxire variables 

a , effective absorptivity of gas for black-body radiation at tem- 
perature 

P Pl«w/(^+^Pl«wAg) 

a rTg^/Bu kg 
Pi a nr^^/Nu kg 

7 ratio of specific heats 

7 ratio of specific heats of air at NACA standard sea-level condi- 

° tions, 1.40 


€ 


effective emlssivity of gas 
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emisslvity of wjxe 
(4 Nu kg/D^ 

6 intermediate time variable 

K thermal diffusivity of wire 

II gas viscosity 

I intermediate distance variable 

p density 

0 Stefan- Boltzmann constant 

T time constant, T^/(i+4"^€^€g/T^) 

Ti 0^4 Nu kg 

Tpt value of for platinum wire 

Tq value of Tpj. at NACA standard sea-level conditions 

$ transient tenqperature function 

y\f spatial temperature distribution function 

to angular frequency 

Subscripts : 

A refers to wire material A 

av arithmetic average 

B refers to wire material B 

b value at wire support 

g refers to gas 

1 value at x =» (L-L')/2 and at x = (L+L')/2 

m value at x ■ l/2 
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P refers to wire diameter Dp 

Pt value for platinum 

Q refers to wire diameter D 

Q 

s refers to static conditions 

t refers to total conditions 

w refers to wire 

0 refers to MCA standard sea-level conditions 

1 initial value 


2 


final value 
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APPENDIX B 

MATHEMAiriCAL SOLOTIONS OF HEAT-TRANSFER 
EQUATION FOR SOME BASIC DESICaJS 


I. General Equation 


Subject to the approximations and assumptions listed under THEQRT, 
the general equation (lO) is derived from a consideration of the balance 
between the rate of heat transfer to an element of length of the wire 
and the rate of storage of heat in the element. The general equation is 
repeated here for convenience; 


T 



1 



(Bl) 


For consideration of the transient response, it will be convenient to 
transform equation (Bl) into 


5t Sx^ 


(B2) 


by a shift of the origin of ordinates, 


y s T^ - Tf 


(B3) 


so that y represents the difference between Instantaneous wire tem- 
perature and the temperature which the wire would ultimately assume in 
the absence of conduction. 

In general, T and y are functions of x and of t. The 
values of these two tenq)eratTares at the midpoint of the wire, at wMch 
the thermocouple junction is located, will be denoted by ^ and yjj^, 

respectively. The characteristic parameters T and q, plus the ini- 
tial and boundary conditions, then serve to define fully the values of 
T^(x,t), y(x,t), T^ jj^(t), and y^Ct). The conduction parameter q 

governs the steady- state temperature distribution; the time constant T 
exercises the principal influence on the transient behavior, although 
the transient behavior is also influenced by q. 
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II. Steady-State Temperatiire Distribution in 

Single Wire Stretched between Supports 

If both thermocouple materials have essentially the same diameter 
and thermal conductivity, the physical situation shown in figure 2(a) 
may be treated as though a single wire, possessing characteristic param- 
eters T and n, were stretched between supports at constant tempera- 
ture T, . In the presence of a constant gas temperature T , the 
steady state of temperature distribution in the wire is then represented 
by setting ST^/St equal to zero in equation (B1), so that the equation 
becomes 

1 Tv , ^ 

^ + T. - T = 0 (B4a) 

n2 dx2 ^ 

subject to the boundary conditions 

T » T, at 3C = 0 (Mb) 

w b ^ ' 

at X = L (Me) 

The solution of this equation is 

Tw = Tf + (T^ - T^)a|/ (B5) 


where 


Vi; = 


sinh Tyx + sinh (ryL - tix) 
sinh tiL 


(B6) 


In particular, the value of 
X = l/ 2 will be denoted by 


ij/ at the midpoint of the wire 
so that 
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vhere 


= sech (tiL/2) 


(B8) 


Figure 4 shows graphs of against i/l for various values of 
the dimensionless parameter tiL. It is to be noted that only in the 
decade 1< t^KIO is markedly different from the values 0 or 1. 


III. Steady-State Temperature Distribution along Thermocouple 
Wire Pair Stretched between Supports 


In the general case where a thermocouple composed of wires having 
different diameters or thermal conductivities is stretched between sup- 
ports at constant temperature T-^j, subscripts A and B will be 
assigned to represent the two wire materials. The situation is illus- 
trated by figure 2(b). To treat this problem, as well as the one in 
the following section, it is convenient to utilize the following: 

Lemma; The solution of the equation (B4a) for steady-state temper- 
atTire distribution along a single wire subject to the boundary condi- 
tions 


'-W -"p 


at X » Xq, 

(B9a) 

at X >■ Xg 

(B9b) 


is 


(T„ - Tf) sinh (tjx - t)Xq) + (Tg - Tf) slnh (t]x_ - nx) 
Tv = Tf + ^ 

sinh (t)x - Tixg) 
a ^ 


(BIO) 


For the problem of two thermocouple wires, A and B, the differ- 
ential equations and their associated boundary conditions are: 


0 < x< 


L 

2 


For wire A 
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1 Tv,a 


Ha' 


2 dx2 


^ + Tf ,A - Tv,A = 0 


^ 


Tv, A = T^ ^ = 0 


Tw,A * ^ ~ 2 


y 




(Bll) 


For wire B 




1 T, 


v,B 


^b' 


.2 dx2 


^ + T. 


f,B " -^w^B 


■^w^B “ 




Tw,B “ Tv^ni ®'*' ^ “ 2 ^ 


J 


(B12) 


It will be noted that continuity of at x = l/ 2 has already been 

postulated by the preceding expressions, and that the xinknown has 

been used as one of the boundary values. In order to determine ^ ^ 

the additional postulate Is made that the net rate of flow of heat * 
across the thermocouple jxmctlon shall be zero. Since the rate of flow 
of heat along wire A into the Junction is 


^,A 


* d T. 


v,A 


dx 


and the rate of flow of heat along wire B from the Junction is 


4 dx 


it is necessary that 
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?. ^ ^W,A -n 2 ^ ^W-B — 

k , D. “ ^ "R ^ ~ 2 

w,A A (ix “ dr 


(B13) 


Equations (Bll), (B12), and (B13) constitute the complete statement 
of the problem. 


By use of the lemma previously stated, appropriate values for 
and Tfl are inserted into equation (BIO) to obtain solutions of equa- 
tions {Bll) and (B12). Inserting the value z = L/2 into these solu- 
tions and then applying the condition stated by equation (B13) give the 
solution for T as 



where 

“a 

“b “ ^,B ^B ^B 


(B15) 


Expression (B14) may be simplified if the steady-state radiation 
errors for the two wire materials are assumed to be sufficiently similar 
that ^ and g may each be replaced by a common value (for 

example/ “ ("^f A assumption, equation 

(B14) becomes 



(B17) 
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and 


m. 


A,B 


_ 1 


m. 


B 


“b 


(BlTb) 


Tlrus an. equivalent value ^|r m nay Be defined to replace tiie quan- 
tity in equation (B7). Alternatively, an equivalent quantity q’L 

may Be defined to replace tBe quantity tjL in equation (B8), where q' 
will lie Between and A taBulatlon of q'L as a function of 

the dimensionless quantities B gl'^’en- in 

table III. ' 

The experiments descriBed in this report show that, for the range 
of Beynolds and Mach nuniber covered in these experiments, the Nusselt 
numher varies directly as the square root of wire diameter at any given 
set of gas-flow conditions. Hence, if radiation is neglected so that 
T] may Be replaced By in equation (llf), the application of equa- 

tions (lie) and (B15) leads to 



1.25 


(B18) 


Values of (k^ are. given in table II for standard thermo- 

couple wire pairs to facilitate computation of m^ j. 

It is evident that if Both wires have substantially the same diam- 
eter and thermal conductivity. Both - q^) and m^ ^ close to 

zero, and \|/'m given By equation (B17a) is substantially the same as 
\1/^ given By equation (B8) . 

Typical temperature distribution patterns along the wire for vari- 
ous values of ^ case whore q^^L = 4 and q^L » 10, are 

shown in figure 20. It is arvident that, depending on the value of m^ g, 
the thermocouple Junction temperature will lie somewhere Between the '' 


values of T, 


w,m 


for each of the two wires as given By figure 4 or 


equations (B7) and (B8) . 
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IV. Steady-State Temperatiire Distribution along Wire Stretched 
between Intermediate Wire Supports of Different Diameter 

If, as shown in figure 2(o), the wires constituting the thermo- 
couple Junction are supported by intermediate wires of different diam- 
eter, the intermediate wires also exposed to the gas stream but with 
their bases at fixed support temperature T^, the values of tj will be 
different for the intermediate wires even though they are of the same 
material as the adjacent wires forming the Junction. 

For simplicity, it will be assumed that both intermediate wires 
have the same diameter, both wires forming the Junction have the same 
diameter, and both material A and material B have the same thermal 
conductivity. (The results of the previous section can serve to correct 
for deviations from these assun 5 >tions . ) Under these assumptions, the 
intermediate wiires, each of length (L - L')/2, may be assigned an t) 
value of Tjp and the two wires forming the Junction may be treated as a 
single wire of length L’ with an tj value of the subscripts P 

and Q now seiving to distinguish between the two diameters. By virtue 
of symmetry, the temperature ^ of the midpoint of wire Q is 
closer to the effective gas temperature Tj q than is the ten^ierature 
of any other point of the assembly. Similarly, for reasons of symmetry, 
the temperature of each Junction between the central wire and the inter- 
mediate wires may be represented by the common symbol Tv, i* 

A line of reasoning similar to that of section III regarding the 
rate of heat flow across the Junction between central and intermediate 
wires leads to the following mathematical statement of the problem: 





+ 



■•w,? 


0 ^ 


Tw,P - Tt at X = 0 \ 


’■w,P • ’■v,! 

, ^ L-L* 

at X^ — 

j 


^f,Q ■ '^w,Q = 

0 

/ 

'^w,Q ” ^w,l 

, L-L* 

at X 

V 

’^w,Q “ “^w,! 

* , L+L* 

a*t X = 2 

J 


0 < X < 


L-L* 

2 


L-L* 

2 


< X < 


LfL* 

2 


(B19) 


(B20) 
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dz 



dx 


at 


z 


L-L« 

2 


(B21) 


Tlie presence of synmetry makes it unnecessary to state or to use the 
relations for x > (L+-L’)/2. They would he similar to equations (B19) 
and (B21) . 


Use of the lemma stated in section III, Insertion of the boundary 
conditions from equations (B19) and (B20) into equation (BIO), and use 
of equation (B21) to eliminate the intermediate quantity ^ yield 
the following solution for ’ 



(B22) 


where 

^ ^ ^P 


(B23) 


Equation (B22) may be simplified if it is assumed that the effects 
of the steady- state radiation errors are sufficiently alike that p 
and Tf q on the right side of the equation may each be replaced by^a 
common value Tj (for example, » (T^ p + T^ q)/ 2). Under this 

assumption, equation (B22) becomes ’ 


T 

w,m 


Tp + 

f,av 


(T> 




»T 

m 


(B24) 
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Thus oquivaXeirt values jjj and t)”L may "be so defined as to 
replace the quantities and t^L in equations (B7) and (B8). A 

tabulation of q"! as a function of the dimensionless quantities 
qp(L-L ’)7 (niQ/®p) is given in table 17, Because both 

materials have been assumed of the same thermal conductivity, equation 
(bis) and the discussion in section III pertaining to this equation show 
that nu/nip • is given very simply by 



(B26) 


By similar reasoning, it can also be shown that 





(B27) 


Temperatttre distribution patterns for a few cases of practical 
interest are shown in figure 21 . 


7. Transient and Steady- State Eesponse of Single Wire to Step Change 
in Gas Temperature, Wire Initially at Support Temperature 

This case, because of its relative mathematical simplicity, pro- 
vides a convenient introduction to the subject of transient response of 
the wire. The general equation (B2) and its corresponding boundary and 
initial conditions for this case are 


^ -,2 3^2 
y * T^ - T^ at x = 0 

y ■ T^ - T j at x = L 

j - T-^ - at t » 0 


y = - ^f 


> (B28) 




Asstnne a solution of the form 


y » u + V 


(B29a) 
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where u = ufi), with 


1 - 

u - - Tj 

u » 

and V = v(x,t), with 



1 

T ~ - 

at “ 

? 

V = 0 

at 

V = 0 

at 


A 


u = 0 




at X =» 0 



(B29b) 

at X = L 

J 



a^v 




X = 0 
X » L 

1 


(B29c) 


V ai T-jj - - u(x) at t 



The solution of equation (^29h) has already been derived in sec- 
tion II. This solution provides the steady-state solution for j, 
which is rewritten here as 


u » (T^ - Tj)>|r (B30) 

where \|/ is defined by equation (B6). 

The solution of equation (B29c) will provide the transient response. 
Solution of the equation is facilitated by the substitution 

w » w(x,t) = ve"*/^ (B31) 


which converts equation (B29 c) into 
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Ti 2 T — 

_ 


at 

ax^ 

V 

» 0 at 

X = 0 

V 

a 0 at 

X = L 

(x) 

II 

Tj)(l - \|;) at t 


The solution of equation (B32) is (reference 9, paragraph 34) 



1 1 
I (VlfKl-M-) d{ e ’ 


Carrying through the integration and then returning to the variable v 
yield the solution of equation (B29c) as 


V = (T^ - Tf) 9 


(B34a) 


where 


$ » — e 

7C 


- t/r 


2 2^ 
n * t 


l-cos(nx) 

n - ._2 




l+rinlnt 


The term - cos (n*^ makes the series alternating, with only odd 
values of n present. 

Addition of the expressions for u and v to obtain y (equation 
(B29a)) and replacement of y by (T^ - Tf) yield the general solution 


"^w " '^f ^"^b " ^f + (^b ■ "^f ^ * 
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where t|/ and $ are given by equations (B6) and (B34b), respectively. 
At the midpoint of the wire, 2 : = L/2 and equation (B35) becomes 


= ''f + <’^6 - > *m 


(B36a) 


where 


(tL/2) 


= - e 

m 


‘ ^ - 

n»l 


^njt'j I 

UnV/n^L^ 


H It ^ 

■ “572“ 

e n L T 


(B8) 


(B36b) 


At t = 0, is equal to (l For t > 0, 4^^ maybe 

represented very closely (to within 0.03) by the approximation 


»a ■ (l-'V) 


(B36c) 


so that an adequate approximation for T, 


w^m 


is 


T ^ 
w^m 


Tf + (T^ - (T^ - Tf)(l 




(B36d) 


The physical significance of this solution has been discussed 
briefly in the section THEOET and illustrated in figure 5. The asymptote 
of the temperature-time curve is Tj. - (Tf. - The initial slope 

is (T^ - T-[j)/t, which is independent of the conauction error factor 
The effective time constant is x(l Thus, when there is appreci- 

able conduction, the indicated magnitude of ultimate temperature change 
is reduced by an amount (T^. - from its "true" value Tf, but the 

time constant is altered in the same proportion in such manner as to 
maintain the initial rate of change equal to that which would be char- 
acteristic of the zero-conduction case. 
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VI. Transient and Steadj-State Besponse of Single Wire to Step Change 

in Gas Temperature, Wire Initially at Gas Temperatiue 

If the wire is initially at a temperature [^f^l + (T-jj - 
corresponding to an actual gas temperature Tg related to T^ by 

an equation of the form of (llg), and a step change in gas temperature 
occurs to a value of Tg^g Tieldlng a corresponding effective gas tem- 
perature Tf^g section II, an ultimate wire ten^terature 

^Tf 2 + (^h " then the mathematical expression of the problem 

calls for solving equation (B2) subject to the boundary and initial 
conditions 

7 = ^b - ^f,2 
y = T-b - Tf^g at 
7 “ Tf,l + (Tb - Tf,i>l' - = Oj 

Assmaing, as in section V, a solution of the form 

y = u + V (B38a) 


X - 0 '> 

X = L \ (B37) 


where 


u = u(x) u»Tj^-T ^2 x = 0 and at x = L 

V = v(x,t) V = 0 at X » 0 and at x = L 

V =. T^^j_ - Tj^2 + (^"b ■ t = 0 

a procedure similar to that of section V yields the solution 
■^w • '^f,2 (’^b ■ ^f,2^'t' + ("^f,! " ^f,2)® 


(B38b) 


(B38c) 
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wMch at X = l/ 2 becomes 

Tv,m - + Cfb - Tf.ajV ♦ (Tf,l - Tf,2)*lE C“°) 

where $, and are defined by equations (B6), (B34b), (B8), 

and (B36b), respectively. 

In accordance with the approximation (B36c), the physical signifi- 
cance of the solution is illustrated by figure 3(a) . The Initial value 
of the indicated temperature-time curve is PfA- (Tb - the 

asymptote is + ('^b " J the indicated magnitude of ulti- 
mate change is (Tf^l - Tf^2)(^ initial slope is 

(Tf ]_ - Tf 2 ) A independent of the conduction error factor \^ui; 

the effective time constant is T (1 - . Thus, when there is appreci- 

able conduction, the indicated magnitude of ultimate temperature change 
is reduced by a factor (l - "but the apparent time constant is 

reduced in the same proportion, so that the initial rate of change 
remains the same as in the zero-conduction case. 


YII. Transient and Steady-State Eesponse of Single Wire 


to Suddenly Applied Sinusoidal Variation in Gas Temperature 

If the wire, initially at temperature + (T-jj - Tf^av)^^ 

subjected, at time t = 0, to a sinusoidal variation T sin (ciJt) in 
gas temperature superimposed on the steady-state value T^ the 

differential equation with its boundary and initial conditions is 


ar„ 1 / ^^T, 


^ St 

*?/ 




T^ » T|j at X =» 0 



Tw = Tb at X = L 

T 

w 


+ (T. - )\|/ j 

av b f,av 


- T 




(B41) 
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Introduction of the transformations 


T = - T sin (ort) 

z = T e -T, __ j — 




(nji^ \ t/rr —1 

e jcos (ci)t) + (OT sin ((ot) 


(B42) 

(B43) 


changes equation (B41) to 


— ■ A. 

St 




sin((ob^ ^ e^^^os(<it)+(OTsin (0)t^ 


(B44) 


at X s 0 and at x ■■ L 


z * (Tv - T, H 


T(i)T 


at t « 0 


J 


The solution of equation (B44) is (reference 9, section 34) 

-L 


r,2^2+ 

n ir u 

2^2 

q L 


\ ^ ^ /‘’’ll - cos (mr) - Tf^a^ - T sin {(iff) + 


nit 

tI^Lt 


COT 




(cos (06 + COT sin C06) 



(B45) 


For convenience. 


write 
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i s ^ [l - cos (nx)1 sin 
n njt j_ J ^ 

B X- + n^ 


Ttien, performance of tiie integrations in equation (B45) and return to 
the variable T-,, by use of equations (B42) and (B43) result in 


Tv ■ 


V .2 2 2t- 2 r ,/ VI 

— \ A +(j) T B -1, , -l/cOTBi — 

T , o o ” t*^-'tan (0JT)+tan \ —r — ) +TwtC 

V 1+of L \ ^ / 


where 


. , ^ MVD 


^ a -5sJli. 


^ "" b 2+<^'r2 


and use has been made of the identities 


) 

Z / b 

n=l n 
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At the midpoint of the. wire^ vhere x * L/2, equation (B48) "becomes 


vhere 


^w,m " ^f,aT m 


■ Vii^WTV 

Vl+ca^T^ 


sin |wT- tan“^ (C0T)+tan“^^^^J|^+T(0TC^ (B51a) 




A . 1 - 

L f ^ h ^4«a‘ 

nsl ^ 

0 » 

a S — : 

h 2 + 0 ^. 

n=al 


J-fi 


V' 

/ ®ti,m 2 2~ 

n?=l 


m 


(BSlh) 


(B51c) 


(B51d) 


^ m " (W2) 


s i. H 

n^m U 3 t 2 


(B8) 


Approximations (to about 0.03 in anqjlitude and to about 3° in angle) 
for the functions appearing in equation (B51) are 


EZ^5S= - V 

\ 1 + ti^T^ + (a^T^(i - 


tan“^£0T - tan”^ 
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- e ’‘'“3/ [n. mV (1 .,1, 


SO that an approximate solution is 




T„,„ - Tf,e,+ (Tb-Tf,av>4'm-'' , ~ .j aln^-ten-^[mT 


/I l - Ki ) 2 T . 2 ( i .^|,^)2 


(lA) 

T tor — ^ — ® 


•V[t(i-'I-„^ 


(B53) 


Thus, hy comparison with aquation (3), the effect of conduction is to 
make the system act as though it possessed a time constant t( 1 
rather than T. At the same time, the amplitude of the steady- state 
tezna and the amplitude of the transient term are hoth attenuated hy a 
factor (l “ . These conclusions are similar to those obtained for 

the case of a step change in effective gas temperature. 
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APPEHDIX C 


HUMERICAL COMPUTATIONS 

The time constant and the steady- state radiation and conduction 
error may be expressed, to an accuracy sufficient for all practical pvir- 
poses, in terms of the experimentally measured static pressure and total 
tenperature and of the Mach number, the value of which is readily 
obtainable from the experimentally measured static and total pressures. 
The numerical formulas for the intermediate quantities required for the 
numerical expression of the approximate values of radiation error, time 
constant, and conduction error will be derived in sequence. 

Viscosity . - The viscosity n of air is independent of pressure 
£ind can be expressed in terms of the tenperature by means of the formu- 
las 


h - ^0 



(Cla) 


H - 0.159X10"® T°’®® lbm/[(sec)(ft)] 


(Clb) 


where T is in °E, to within 1 percent in the range 460° <T <2300° R, 
if values given by reference 3, quoting a personal communication from 
Genereaux, are assumed. In equation (Cla), Uq is the viscosity of 
air at temperature Tq. Equation (Clb) is based on assuming equal 

to 11.9X10“® pound mass 
be assumed that formula 
400°< T<3000° R. 


per second per foot at Tq « 519° R. I? will 
(Clb) can be extrapolated to the interval 


Reynolds number . - The Reynolds number Re can be expressed in 
terms of Ifech number, pressure, temperature, and wire diameter by use 
of formula (Cl) for the viscosity, the universal gas law, and the for- 
mula for the velocity of so\md 


where r 
written 


Re 






is the specific gas constant for air. 


D ^ (C2a) 

Formula (C2a) may be 


Re » l.OlxlO® DMpTg"^-^ 


(C2b) 
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where D is in inches, p is in atmospheres, and T is in °R. The 
Reynolds number so defined involves the evalmtion of the gas properties 
at static pressure and static temperature. If, in order to improve the 
correlation between heat transfer and aerodynamic conditions, a differ- 
ent Reynolds number Re* ,ia defined that Involves evaluation of the gas 
density and viscosity at total rather than at static temperature, the 
conversion formula 

^'t “ ’^s ^ (C3) 

may be applied to equation (C2) to yield 

Re* - — ( -2] M i/ytT 7 D -2- » Re / (C4a) 

^io\V ® ^t 

This formula may be written 

r / 

Re* - 1.01X10^ D M p ^ (1 + j (C4b) 

if D, p, and T are in inches, standard atmospheres, and °R, 
respectively, and Xq (•■ 1*40) is the value of y at NACA standard 
sea-level conditions. 

Thermal conductivity . - Because adequate experimental data on the 
thermal conductivity of air at hi^ temperatures are unavailable, the 
thermal conductivity kg will be computed on the assumptions that the 
Prandtl number is constant and equal to 0.71, the viscosity is given by 
equation (Cl) , the specific heat of air is given by the values tabiilated 
in reference 10 for the temperature range 460°<T<4000° R, and the 
specific heat - tenq)erature relation can be extrapolated to 400° R. 

Since little Information is available regarding th^ variation of Prandtl 
number of air with temperature, and some of the information is contra- 
dictory, it has been recommended by Dr. E. R. G. Eckert of the Lewis 
laboratory that an assumption of constant Prandtl ntmiber be made. The 
value suggested by Dr. Eckert is the latest value at room temperature, 
as given by reference 11. The formulas 

kg - kg^o (C5a) 

kg - 3.03X10"® T°*^® Btu/fft)(sec)(OR^ (C5b) 
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where T is in °E, then fit the computed values within 2 percent in the 
range 400° <T< 3000° R. In formula (C5a), hg^O is the thermal con- 
ductivity of air at temperature Tn- In form^a (C5b) , the value of 
kg,o is taken as 4 .0X10"® Btu/[(ft) (sec) (°R)J at Tq » 519° E. 

Nusselt number. - Insertion of formulas (C4a) and (C4h) into equa- 
tion (26a) lesids to 


Hu 


0.431 



M p 



(C6a) 


Hu 


1/2 -0.60 

1.37X10^ (D M p) 



(C6b) 


where, in equation (C6h) , the units of D, p, and are inches, 

atmospheres, and °R, respectively. 


The computations of the quantities T^, and t)]_ require 

knowledge of the thermal conductivity of air. The uncertainty in knowl- 
edge of this quantity is sufficiently large that no appreciable increase 
in total error resiilts by xise of y « 1.34, as an average value over 
the temperature range 500° < T-t < 3000° R, for computation of 
anfl Tjj^. Formula (C6b) then becomes 


Nu 


1.36X10^ (D M -r (1 + 0.2M^) 


(C6c) 


Steady-state radiation error . - Computation of the steady-state 
radiation error is conveniently accomplished by first evaluating the 
quantity . Insertion of equation (C5) (evaluated at total tempera- 
ture) and equation (C6) into equation (l3d) and replacement 

quantity by the almost numerically equal quantity Ty yield 



2.32 a 





(C7a) 


10^ 3.82, 


Pi - 0.97X10 * Ty 



+ 0.2m2) 


1/4 


(°R) 


(C7b) 
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where, in equation (C7b) , the units, D, p, and are inches, 

atmospheres, and °R, respectively. The nomograph of figure 16(a) yields 
when M, D, p, and 1^. are given. 


The radiation error may then he computed as 


where 


Radiation error = 









1 + ^Pl^w"g/T^ 


Pl^w 


(13b) 


(13c) 


Tn the absence of conduction, the actual gas temperature Tg is then 
obtainad by subtracting the radiation error from the indicated wire 
temperature Ty. 

The quantity represents approximately the error present when 

T^/T^r ■ 1.2. As an aid to confutation, a nomograph for (T^/T^)^ is 
presented in figure 16(b) . 

constant . - The computation of time constant in the 

absence of radiation or conduction, is accaif llshed by inserting equa- 
tions (C5) and (C6) into equation (lla), yielding 


1/2 1/4 0.435 

0.5S1.0 ' r ' Tp _ 

'^ 8,0 




3/2 


. ,1/2 0.18 1/4 

(M p) ' T^ r ' 


1/4 

1 + ^M^) (C8a) 


3/2 


,1/2 0.18 


p c D ^ 1/4 


(1 + 0.2M2) ^ sec 


(C8b) 


(M p)^' T^ 


where ^ in equation (C8b ) , the units are 
p^c^ Btu/ [(ft^)(°H)] 

D inches 

p atmospheres 
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A nomograph for the calcxilation of appears in figure 17. The 

first turning line has been graduated to give directly the value 
the time constant of a platinum vire at NACA standaurd sea- level condi- 
tions. The final result of the nomograph is the time constant 

of a platinum vire at the given aerodynamic conditions. The time con- 
stant '’’i for vire of any other material is obtained by multiplying 
Tpt by (p A Ppt^Pt^ ■ constant T in the presence of 

radiation is 


T - T^/(l + 4Pi€^/Tg) « Ti/(1 + 4^Cv€g/T^) (C9) 

A tabulation of (p^c^/Pp^Cp^) for various thermocouple materials 
is given in table I. A tabulation of the arithmetically averaged values 
for commonly used thermocouple pairs is given in table II and is 
repeated for convenience in figure 17. The values appearing in table I 
vere obtained from references 12 and 13 and are believed to represent 
the best values available. The data in figure 17 are thus sxifficlent 
for estimating 7^ for any group of aerodynamic conditions that may 

normally be encountered in aircraft engine operation. 

Steady-state conduction error . - Computation of the steady-state 
conduction error is conveniently accomplished by first evaluating the 
conduction parameter nL. If the time constant T has already been 
determined, t]L may be computed most simply by use of equation (ill) as 

(=10) 


The diffusivities of various materials are listed in table I. 

If the time constant is not knovn, t]L may be computed by insert- 
ing equations (C5) and (C6) into equation (He), yielding 


1.72 k 


(^lL)‘ 


g.Q 

1/2 1/4 0.435 


(M p) 


1/2 0.18 
' m r 


3/2 
D ' k 


w 


t1/4 




P -3 (D M p)^/\°*^®(L/D)' 

(tj^^L) » 1.65X10 ^ 


kv (1 + 0.2M^)^/^ 


(CHa) 


(CHb) 
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where in equation (CUb) the units are 
Btu/ j^aec)(ft)(°E^ 

D,L inches 


p atmospheres 

A nomograph to facilitate calculation of presented in fig- 

ure 18. The final result of the nomograph is (q;,L)p+ , the value of 

2 -l- A U 

( H2_L) for a platinum wire at the given aerodynamic conditions . The 
value of for a wire of different material is obtained by 

multiplying , Values of (kj^/k^) are presented 

in table I. A tabiolation of the harmonically averaged values of 
for commonly used tb.ermocouple pairs is given in table II and is 
repeated for convenience in figure 18. The value of (tL)^ is then 


(1 .• 4Bi<yTg) » (1 1 C 012 ) 


Finally, the conduction correction factor is 

. sech (qL/2) (B8) 

and is given graphically in figure 19. 

Corrections for differences between thermocouple vires . - If the 
two thermocouple materials differ appreciably in thermal conductivity, 
the computation of (^i^) formula (Cll) should be carried out for 
each half of a wire, and then the effective for the pair of wires 

deteimined as indicated in part III of appendix B . If the two thermo- 
couple materials are sufficiently alike, the effective thermal conduc- 
tivity is the harmonic mean of the values for the individual materials. 

If the two thermocouple materials differ in the values of and 

of "til® effective values of these quantities are the arithmetic 

means of the respective values for the individual materials . 

constant in presence of radiation and conduction . - The time 
constant in the presence of radiation and conduction is 
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1 

T (1 > i|f ) « T — S — 

“ ^ (1 + 4P^e^Cg/T^) 


(C12) 


Gas temperature In presence of radiation and conduction eiTors . - 
The gas temperature when the thermocouple is subject to both radiation 
and conduction errors is 


g w 


/m \4 

(l-a J (^) - (1-eJ 

^ \Tv/ S 




r\ 


T. - e - (l-£g) - (VV ♦m 


(C13) 


where Tj^ is the support teniperature . 
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TABLE I - .PROPERTIES OP THERMOCOUPLE ELEMENTS AND ALLOYS 



Platinuni 1334 0.0324 0.01142 0.0381 1 1 

m^odiuni 774 . 058 . 0125 . 040 1.05 .92 

Platinum 1261 .0357 ^ .OO 484 -0155 1.05 2.38 

plus 
13 per- 
cent 
rhodl'um 

ALumel 537 .124 ® .0048 *0103 2.3 

Chromel 545 .106 ® .0031 .0077 1.34 3.7 

Constantaua 553 . 099 . 0038 .OlO^^ 1.27 2.9 

Iron 491 .107 .0096 .026 q 1.22 1.19 

Copper 555 .093 .O 6 I 0 .172 1.20 .19 

Aluminum 169 .220 .032s «126 0.87 .35 

^Evaluated for range 530° to 660° R except vhere noted. 

Evaluated at 530° R. 

°Evaluated at 660° R. 
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TABLE II - PROPERTIES OF THERMOCOUPLE PAIRS 



[subscript A denotes positive material. This material is listed 

first in the palr^ 


Material 

CPw°w)A-^CPw=v)a 

^Pi / 1 1 r \ 


■Mns 



^ V\,A \,W 

(vJ 

imiii 

Platinum plus 
13 percent 
rhodium - 
platinum 

1.02 

1.69 

1.54 

0.0220 

Chromel-alumel 

1.45 

3.0 

.81 

.0088 

Chromel- constantan 

1.31 

3.3 

.90 

.0087 

Iron-constantan 

1.25 

2.1 

1.59 

.0146 
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TABLE III - VALUES GF tj'L 



tj’L 

^ A ^ 

1 

"" a^b 


1 

2 

3 

4 

5 

10 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.00 

1.00 

1.00 

1.00 

2.00 

1.66 

1.40 

1.19 

1.00 

3.00 

2.16 

1.67 

1.30 

1.00 

4.00 
2.54 
1.85 
1.30 

1.00 

5.00 

2.80 

1.97 

1,43 

1.00 

10.00 

3.25 

2.15 

1.51 

1.00 

2 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.19 

1.40 

1.66 

2.00 

2.00 

2.00 

2.00 

2.00 

2.00 

3.00 
2.69 
2.43 
2.21 

2.00 

4.00 
3.25 
2.73 
2.33 

2.00 

5.00 
3.68 
2.94 
2.41 

2.00 

10.00 

4.52 

3.28 

2.54 

2.00 

3 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.30 

1.67 

2.16 

3.00 

2.00 

2.20 

2.43 

2.69 

3.00 

3.00 

3.00 

3.00 

3.00 

3.00 

4.00 
3.70 
3.44 
3.21 

3.00 

5.00 

4.27 

3.75 

3.34 

3.00 

10.00 

5.55 
4.30 

3.55 
3.00 

4 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.38 

1.05 

2.54 

4.00 

2,00 

2.33 

2.73 

3.25 

4.00 

3.00 
3.21 
3.44 
3.70 

4.00 

4.00 

4.00 

4.00 

4.00 

4.00 

5.00 
4.70 
4.44 
4.21 

4.00 

10.00 

6.48 

5.28 

4.54 

4.00 

5 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.43 

1.97 

2.80 

5.00 

2.00 

2.41 

2.94 

3.68 

5.00 

3.00 
3.34 
3.75 
4.27 

5.00 

4.00 
4.21 
4.44 
4.70 

5.00 

5.00 

5.00 

5.00 

5.00 

5.00 

10.00 

7.33 

6.23 

5.52 

5.00 

10 

1.0 

.5 

0 

-.5 

- 1.0 

1.00 

1.51 

2.15 

3.25 

10.00 

2.00 

2.55 

3.28 

4.52 

10.00 

3.00 

3.55 
4.30 

5.55 
10.00 

4.00 

4.54 

5.28 

6.48 

10.00 

5.00 

5.52 

6.23 

7.33 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 
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TABLE IV - VALUES OF n"L 
_ 


mp 

Tjp ( L - L ') 

“Q 

1 

2 

3 

4 

5 

10 

0.01 

1.45 

2.31 

3.27 

4.26 

5.25 

10.2 

.02 

1.45 

2.32 

3.28 

4.27 

5.26 

10.2 

.05 

1.47 

2.34 

3.30 

4.29 

5.29 

10.2 

.10 

1.51 

2.38 

3.35 

4.34 

5.33 

10.3 

.20 

1.57 

2.46 

3.43 

4.42 

5.42 

10.4 

.50 

1.75 

2.68 

3.66 

4.66 

5.66 

10.6 

1.0 

2.00 

3.00 

4.00 

5.00 

6.00 

11.0 

.01 

2.31 

3.04 

3.94 

4.90 

5.89 

10.8 

.02 

2.32 

3.05 

3.95 

4.92 

5.90 

10.9 

.05 

2.34 

3.09 

3.99 

4.96 

5.94 

10.9 

.10 

2.38 

3.15 

4.06 

5.03 

6.02 

11.0 

.20 

2.46 

3.26 

4.19 

5.17 

6.15 

11.1 

.50 

2.68 

3.57 

4.53 

5.52 

6.51 

11.5 

1.0 

3.00 

4.00 

5.00 

6.00 

7.00 

12.0 

.01 

3.27 

3.94 

4.80 

5.76 

6.74 

11.7 

.02 

3.28 

3.95 

4.82 

5.77 

6.75 

11.7 

.05 

3.30 

3.99 

4.87 

5.82 

6.81 

11.8 

.10 

3.35 

4.06 

4.95 

5.91 

6.89 

11.9 

.20 

3.43 

4.19 

5.09 

6.06 

7.05 

12.0 

.50 

3.66 

4.53 

5.48 

6.47 

7.46 

12.4 

1.0 

4.00 

5.00 

6.00 

7.00 

8.00 

13.0 

.01 

4.26 

4.90 

5.76 

6.70 

7.68 

12.6 

.02 

4.27 

4.92 

5.77 

6.72 

7.70 

12.7 

.05 

4.29 

4.96 

5.82 

6.78 

7.75 

12.7 

.10 

4.34 

5.03 

5.91 

6.86 

7.84 

12.8 

. 20 - 

4.42 

5.17 

6.06 

7.03 

8.01 

13.0 

.50 

1 4.66 

5.52 

6.47 

7.45 

8.44 

13.4 

1.0 

5.00 

6.00 

7.00 

8.00 

9.00 

14.0 

.01 

5.25 

5.89 

6.74 

7.68 

8.66 

13.6 

.02 

5.26 

5,90 

6.75 

7.70 

8.67 

13.6 

.05 

5.28 

5.94 

6.81 

7.75 

8.73 

13.7 

.10 

5.33 

6.02 

6.89 

7.84 

8.82 

13.8 

.20 

5.42 

6.15 

7.05 

8.01 

8.99 

14.0 

.50 

5.65 

6.51 

7.46 

8.44 

9.43 

14.4 

1.0 

6.00 

7.00 

8.00 

9.00 

10.00 

15.0 

.01 

10.2 

10.8 

11.7 

12.6 

13.6 

18.6 

.02 

10.2 

10.9 

11.7 

12.7 

13.6 

18.6 

.05 

10.2 

10.9 

11.8 

12.7 

13.7 

18.7 

.10 

10.3 

11.0 

11.9 

12.8 

13.8 

18.8 

.20 

10.4 

11.1 

12.0 

13.0 

14.0 

18.9 

.50 

10.6 

11.5 

12.4 

13.4 

14.4 

19.4 

1.0 

11.0 

12.0 

13.0 

14.0 

15.0 

20.0 
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W 'jsqmTO qo'wiii OT8J^.B-a3Ja 


Figui'e 1. - Comparison of Mach and Reynolds number ranges 
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Air Flow 



(a) Wires with same D and tj. 
Air flow 



(b) Wires with different D and tj. 



(c) Wires supported between Intermediate 
diameter struts . 


Figure 2. - Basis thermocouple designs (deyeloped views). 



* Initial radiation error 
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and themal conduotlTltj, 
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Figure i. - TeoQ>erature dlatrlbutlon along alaiple vlre in steady state. 






Temperature 
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'■ ' I Final gae temperature, Tg^g 

6[(l-ag ^)(Tcl/Tw)^ - (l-€g)] = radiation error 

J. i_. 

ta * Final conduction error 

I Final thermocouple temperature 



— — — — — — — — — Duct vail temperature, 



Time, t 

Figure 5. - Reaponse to etep change from support temperature. Thermocouple 

initially at support temperature . 



Temperature 
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Initial tauperature 
due to beating 


|l[(l-ag^(i)(Td/Tv)^ - (l-€g)^ » Initial radiation 

1 " 


” Initial conduction error 
.Jl — — Initial thermocouple temi 


Initial thermocouple temperature, 


-Thermocouple 

response 


■ ■ Gas temperature, Tg 

{?[(l-Og,d)(Td/Tv)* - “ Final radiation error 

— _ — _l_ — — '^f^2 

— = Final conduction error 

Final thermocouple temperature 


Support temperature, 


Duct vail temperature, 


Time, t 

Figure 6. - Reaponae to step change with wire Initially heated to temperature 

greater than gas temperature. 
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Mercury manoTaeter 
(total pressure) 


plenum 

chamber 


-Mercury manometer 
(velocity pressure) 

^Working region 


^ Receiver 


“Vres^sectloT^ 

Thermocouple 


Total- 

temperature 

indicator 


Double pole, 
double throw 
snap- action switch 


Variable 

transformer 


Recording 

oscillograph 


Figure 7* - Setup an d instrumentation for determination or time 
response of thermocouples. 
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Configuration I 


Configuration II 


Metal 


n 


Ceramic 


Taermocouple 

Junctlou 


Probe Thermocouple material 


C hr omel -c ons tantan 0.02 0.15 0.05 0.0085 


.30 .10 .0176 


Iroti-c one tantan 


C hr omel -edumel 


Platinum plus 13- 
percent rhodium - 
platinum 


C hr omel -c ons tantan 


Configuration 


.50 .25 .0554 


.30 .10 


.30 .10 


.0174 


.0176 Straight wire 


^Dimension d is the average of eight measurements . 

Figure 8. - Thermocouples used for test purposes , Air flow perpendic- 
ular to plane of thermocouple loops. (All dimensions are in inches.) 
















































Deflection of recorder. In 
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(a) Chromel-conatantan thermocouple. 


Figure 10. - Time constants of thermocouples with Reynolds number based 
on static temperature. Lines drawn through data points for M<0.5. 
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(d) Pro^je Dj iron- const antanj diameter. 



(e) Probe E; chromel-alumeli diameter, 
0.0155 inch. 



Free- stream Mach miniber, M 


(f) Probe Fj platinum plus 13*percent rhodium - 
platinum; diameter, 0.0174 inch. 


Figure 11. - Concluded. Time constants of seven 
thermocouples . 





Time constant, sec 
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Figure 12 . - Time constants of thermocouples with Reynolds number based 

on total temperature. 









Figure 13. ~ Correlation of Nuaaelt number with Reynolda number for air for all data. Conatant 

Frandtl number. 













I < 1 1 1 1 1 1 1 1 H l-H-H-H-t- 

K-4 ■ 1 1 1 1 1 ■ ii I i f i nH ■ t ■ ■ M 
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Mach Vlre 



Given:' • 500* R Troeedure: Connect H to D and find Intersection A. Answer: • 0.46* R 

p « 1 ata Connect A to Ty> and find Intersection B. 

M - 0.3 Connect B to p and find answer 

B - 0.015 In. 

(a) Noaograph for ecnputlng 

Figure 16. - homographs for determining radiation error. 

(A 10- hy 12- In. print of this fig. is attached.) 
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E<tul7alent 

duct 

tenperature 

(^) 

^3000 

iLasoo 

-^2000 
: 11750 

--1500 

-ll250 

-ijliooo 

-j|« 900 
000 
-L 700 

-j- 600 

-L 500 
i 400 


(Ti/Tv)^ 



(b) Nomograph for ccmputlng (T^/T^)*. 

Figure 16, - Concluded. Nonographa for determining radiation error. 
(A 9- by 10-ln. print of this fig. is attached.) 


Thermocouple 

temperature 


400- 


500 


600- 


1250- 

1500* 

1750- 

2000 - 

2500- 


3000- 
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tUab 

QUSbCl* 

M 

0.1 -r 


static 

pressxzre 

P 

(atm) 

0.1-T- 


0.2’ 


o.s- 

0.2 -h 

0.5-1- 
OJ 

o. g 

1^0- 

0.3+ 4 

2’ 

3- 

4- 

5’ 

6. 

7’ 


0.4^ 

0.5 

0.6 


0.7 

0.8 
0.9 
1.0 


i 10^ 

15 
20 


Thermocouple 

material 


Chrcoel- 

coostantan 

1.30 

Xron- 

constantan 

1.24 

Chrooel- 

alumel 

1.44 

Platinum plus 
IS -percent 1 

rhodium - 
platinum 

1.02 


^At rooB tenperature. 


Total 

temper- 

ature 

Tt 

(®R) 

T- 400 
.500 
800-.^, 

¥ 800^ 
““iOOO 

4-1500 
.2000 
-2500 
-1-3000 


"0 
(•ec) 
5 ■ 

4 -jr 

5 JL 

2 -:r 

1 2 ^ 

0.8 

0.6 -f r 

0.5 -i - 
0.4 -ik 

0.3 

0.2 -±- 








0.I3 




. ^ ^ 


0.004- 
0.003-1- 


Vlre 

llameter 

D 

(in.) 
0.1.^ 
0.09-+- 
0.08’+* 
0.07-^ 
0.06-^ 
0.05-1- 

0.04< 

0.03- 


0 . 02 - 


0.1 

0.08 

0.06 
- ^0.05 
._0.04' 

0.03 

0.02 

0.01 

0.006 

0.006 
0.005 
0.004 

0.003 

0.002 0.001 


0 . 01 - 
0.009' 
0.008 
0.007. 
^0.006 • 
0.005 ' 
0.004 ' 

0.003 ^ 


^ 0.002 i- 


Glv’en; 


M • 0.5 
D - 0.006 in. 
p « 1 atm 
- 500° S 
C^^GBel-alumel 


Procedure: Connect M to D and find 

intersection t q 
C onnect Tq to p and find 
intersection A . 

Connect A to and find 
intersection 7^, 

Multiply Tp^ by number in 
table to find 


Answer : t 


Pt 


0.04 sec 


^Pt 


IfiUii 






■^Pt Pt 
0.057 sec 


Figure 17. - Nomograph for computing 
(A 11- by 12 -in. print of this fig. is attached.) 






KACA TN 2599 


79 


Mich 

mnher 


(Vd) 




3t&tlc 
priisurv, 7 
(hta) 



Vln 

dluwtir, 0 

0.09 « 
0.08 -- 
0.07 -1 


aiTan; N m 0.3 ?Toe*div«: Coouict K to D ftod find latcMictloB A. 

5 - 0.003 ia. Coaaict A to L/D and find Intcrsactlon B- 

t/0 • 30 CoaiMct B to p aad flzjd intorsoetloa C. 

p - I ata CooBoet C to 7. and find lotoraoetloa 

• 5D0* B ' ' ‘I 'Pt ■ 

Chroatl-aluati Italttpijr ty v»lu« In tafala to find . 

MCtita IB. - BoaocTi^ fbr eo^uein< 

(A by 13- la. print of thin n«. u attachod.) 


(nil»)pt^ - X 

(n^L)2 - so 
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0 .2 .4 .6 .8 1.0 


Ratio of dlatanca along vire to total length, x/L 

Figure 20. - Temperature distribution along thermocouple consisting 
of two wires of different thermal conductivity or diameter. 



Figure 21. - Temperature distribution along wire with intermediate supports 
as shown in figure 2(c). One-half syTnmatrical distribution is shown. 

L, 0.5 inch? V, 0.1 inch; np(L-L'), 1.4. 


NACA - Ungley Flald, V*. 






A 




